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INTRODUCTION 

Recently it has become clear that the phenomena in the ionospheric P— region 
are determined to a great extent by dynamical processes in the strato-mesosphere 
and lower thermosphere. It is these processes that significantly influence the 
distribution of the minor neutral constituents and the whole ionic composition 
of the middle atmosphere, regulating the ionisation-neutralization cycle of the 
D-region. In this respect much attention is paid to the study of the winter 
anomaly (WA) phenomenon on medium and short radiowaves, in which the 
meteorological character of the lower ionosphere is most prominent. Significant 
experimental data about the variations of the electron concentration, N, ion 
composi t ion, temperature and dynamic regime during WA permit a better under- 
standing of the character of the physical processes in the middle atmosphere. 

WA is recorded as an abrupt increase of absorption L of the middle and short 
radiowaver on different winter days caused by a significant growth of K 
predominantly in the 75-90 km region with a steep gradient in the 80-85 kn 

interval , WA is usually observed at middle latitudes of both the Northern and 

Southern Hemispheres, and is best recorded in the frequency range f. * f cosi * 

1-2 MHz, where is the equivalent frequency of the incident radiowave with a 
frequency f at angle i. This is confirmed by Figure 1 where we have presented 
the results of the absorption measurements in the Sofia ionospheric observatory, 
along the paths 

Fristira-Sof ia (1412 kHz/170 km; f. - 1.1 MHz) 

Greecc-Sofia (4050 kHz/320 km; f) m 2 MHz) 

All measurements refer to a solar zenith angle \ - 78.5° (cos v * 0,2), 

The seasonal variation shown in Figure 1 allows a quantitative estimation 
of WA magnitude by calculating the excess of the winter absorption value, L , 
over the summer one, L g . in this way the L^,/!.,. ratio, obtained in W 

experimental and theoretical ways, will serve as a criterion for comparing 
the different physical models of the middle atmosphere. 

‘ the purpose of the present report is to evaluate the influence of the 
neutral wind on the seasonal variation of the electron concentration N for the 
altitude interval 90 < z < 120 km, where the ratio v. /. . , 0 f the ion- 
neutral collision frequency, and the ion f.yrof req2enc-y , j decreases 
from 40 to 1. CIRA-72 is used as a model of the zonal wind. 

THEORETICAL MODEL 


The distribution of the electron concentration N in ionospheric regions D- 
and E- is described by the balance equation: 


•N 


iN" - div(N v ) 
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where q is the electron production rate, » the recombination coefficient, and 
% the mean electron velocity. Due to the neutrality of electrical charge the 
following equation holds: 
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Figure la. 


Figure lb. 


i.e. ion velocity should be investigated instead of electron velocity. Then, if 
we assume the three-component model of the ionosphere (neutrals, heavy positive 
ions and electrons) it becomes necessary to determine the ion drift velocity as 
a function of the neutral wind, geomagnetic field and r. - v- / If 
seasonal variations of the neutral wind for different altitudes are known and if 
the preliminary model for q and a is determined and conditioned by iono- 
spheric chemistry in different months, then it is easy to find the seasonal 
distribution of the electron concentration. 


If we neglect the pressure gradient and electric field forces, the macro- 
scopic notion of heavy positive ions in a collisional ionosphere with 
geomagnetic field is described by 


dv./dt * v^ n (v - v^) ♦ (e/n^Mv^ x tf] (3) 

where c and n . are the charge and mass^of the ions, the ion drift velo- 
city, v the neutral gas velocity, and B the local geomagnetic field. Under 
quasi-steady conditions the ion velocity is determined from (3) as follows: 
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where R - 3„C and r • - 1 . Fron (4) it follows that V; is strongly 

influenced not only bv the different neutral winds at different altitudes, but 
also by the strong altitude dependence of r.. When r^ >> 1 there is a full 
drag of iens by the neutrals, and v. v. This holds between 60 and 90 km. 

In the height range from 90-100 to lAO km, is of the order of 1, and this 
region is known as the transition region (AXFORD, 1963). There 


(1 ♦ rj 2 ) -1 Cr. 2 v ♦ .r.[v x Cl). 


(5) 


Finally, when r^ << 1, then v ^ ■ (v * >)>. 

The quantity of the ion flow is measured by the convergence of the ion 
velocity. If v^ is given by (5), then 

div v\ « (1 ♦ rj 2 )" 1 (r. 2 div ♦ r^^rot v)) + 

+ (1 + r^ 2 )” 2 (2r -v + (1 - r^Mv x Cl) Vr. ((>) 
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Experiments indicate that the horizontal neutral wind components usually are an 
order of magnitude or so greater than the vertical component, while the vertical 
scales appear to be the sane for the three components. Hence we can take 
v - (v x# v , 0). We use the standard coordinate system in which the x-axis 
points to geomagnetic east, the y-axis to geomagnetic north, and the z-axis 
vertically upward. A horizontally stratified ionosphere is assumed. For 
altitudes 90-100 km, from (5) and (6) follows 
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•V / •* - It /a ♦ r 2 )} , v /,z ♦ [(1 - r 4 -)/(l ♦ r.-) 2 l , v .r / 
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As we are interested in the seasonal course of N for the respective altitudes, 
we have to solve the quas i-stationary equation (1): 

2 

o 1 "* r . 

•'N“ + [ ly /(l «■ r.^)] r. vjt *' — r./z ♦ r.pv^ H - q - 0 (S) 


where the altitude profile of N is taken in the form N(z) * N exp(pz) 
( VEL1N0V et al., 1974). 0 


In the height range 115-120 kn, the full ion velocity or the equation (5) 
is used. Then the redistribution of electron concentration is influenced not 
only by the zonal but also by the meridional wind. Then: 


v xz ■ V (I + r i )1(r i v x + Vy } 
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From (9) it follows that the influence of the meridional wind increases with 
increasing altitude. 


NUMERICAL RESULTS 

We are interested in the influence of the neutral wind on the seasonal 
course of absorption of short radiowaves, where the WA is clearly felt. Our 

attention is drawn by the courses, shown on Figure 1. As some of the radiowaves 

are reflected below 100 km, other between 100 kn and 110 kn, it becomes 
necessary to consider the altitude range from 90 km to 110 km, and to determine 
mainly the influence of the zonal winds on the seasonal course of N. Assuming 
that L ' N, we can compare theoretical and e:,erimental results. Froo C1RA-72 
we obtain a monthly value of the zonal wind and wind shear for altitudes 90, 95, 
100 and 110 kn, for * 40°N. We include them in equation (6) and solve it. 

For altitudes 100 km and 110 km we used data for the meridional wind from HANSON 
ct al. (1961), for January, February and March only. Their influence on N is 
slight. As the absorption is an integral characteristic, in order to be able to 
compare the seasonal course of L.,. ? with the course of N it is necessary to 
add the seasonal courses of N, intluenced by the zonal wind lor the altitudes 
90, 95 and 100 kn (Figure 2). The influence of the zon^l wind on WA is obvious, 

conditioned by the model of i and q. If the latter increases the winter values 

of N over the summer values by 16*, then if we add the wind this increase is 
augmented up to 25*. The experimental results show an increase by 30*. 
Theoretically obtained seasonal course of N in Figure 2 clearly shows the 
oxper imcntnl ly obtained secondary maxinun in Hay and minimum in July, as veil as 
the sharp increases of N from October to November, known as "October-cf feet M . 
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Figure 3. 


The same is done for the second radio path, and influence of N 

by the altitude 110 km is added. Figure 3 shows the seasonal course of the 
average N. The increases of the winter values are about 382, but considerably 
lower than experimentally observed. This may be due to the above mentioned 
model of q and a. The neglect of the influence of the electric field makes 
itself felt stronger at these altitudes. Figure 3 shows well the simmer 
maximum, which is clearly seen on Figure lb. 

The theoretical results compared with experimental results arc encouraging 
and clearly show the influence of the neutral atmosphere dynamics on the 
seasonal variation of the electron concentration, S. 
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